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In a frequency range where a microwave resonator simulates a chaotic quantum billiard, we have 
measured moduli and phases of reflection and transmission amplitudes in the regimes of both isolated 
and of weakly overlapping resonances and for resonators with and without time-reversal invariance. 
Statistical measures for S'-matrix fluctuations were determined from the data and compared with 
extant and/or newly derived theoretical results obtained from the random-matrix approach to quan- 
tum chaotic scattering. The latter contained a small number of fit parameters. The large data sets 
taken made it possible to test the theoretical expressions with unprecedented accuracy. The theory 
is confirmed by both, a goodness-of-fit-test and the agreement of predicted values for those statistical 
measures that were not used for the fits, with the data. 

PACS numbers: 05.45.Mt,24.60.Ky,11.30.Er,85.70.Ge 



Introduction 

Microwave resonators, also known as "microwave bil- 
liards", are ideal systems to study properties of chaotic 
quantum systems [l|-[5| . Most studies have focused on the 
statistical properties of eigenvalues and eigenfunctions, 
especially on tests of the Bohigas-Giannoni-Schmit (sHU 
conjecture. According to that conjecture, the spectral 
fluctuation properties of quantum systems with chaotic 
classical dynamics coincide with those of random-matrix 
ensembles belonging to the same symmetry class. That 
statement holds up to level spacings determined by the 
period of the shortest periodic orbit of the system. We 
comment on that point below. Data are taken by cou- 
pling the resonators via one or several antennas to sources 
or sinks of a microwave power supply. Because of this 
arrangement, microwave resonators can also be viewed 
as open quantum systems, and measurements of the 
reflected and transmitted intensity amplitudes provide 
generic information on chaotic quantum scattering, each 
antenna acting as a single scattering channel (loj . 

In this paper we report on measurements of the com- 
plex transmission and reflection amplitudes of chaotic 
microwave billiards, and on the theoretical analysis of 
such data. For the latter we use the generic approach 
to chaotic quantum scattering based on random-matrix 
theory (RMT). We have used two types of microwave 
resonators. In the first one, time-reversal (7~) invariance 
holds and in the second one, it is violated by placing a 
magnetized ferrite within the cavity. The large set of 
scattering data taken with either device (considerably 
larger than data sets collected, for instance, in nuclear 
physics) allows us to test the RMT approach to chaotic 
scattering with unprecedented accuracy, both for systems 
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that are T invariant and for those that are not. Most ex- 
perimental investigations of chaotic scattering have been 
restricted to measure cross sections rather than individ- 
ual elements of the scattering matrix S. In our setup 
we use a vector network analyzer to actually measure 
modulus and phase of the reflected and of the trans- 
mitted amplitudes and, thus, of individual S'-matrix ele- 
ments. The additional information garnered in this way 
increases the significance of our tests. In addition to test- 
ing the RMT approach to chaotic scattering, we propose 
and test a method for the determination of the strength 
of T violation from data on S'-matrix correlation func- 
tions. This is of particular significance for those chaotic 
quantum scattering systems for which the relevant pa- 
rameters cannot be determined easily by dynamical cal- 
culations such as the semiclassical approximation [Tl| . 
Some of our results have already briefly been reported in 
Refs. [3G3. 

For T-invariant systems statistical cross-section fluc- 
tuations have been thoroughly investigated experimen- 
tally and compared with theoretical predictions in the 
regime of isolated nuclear resonances [14j | (average res- 
onance spacing d very large compared to average reso- 
nance width r) and in the Ericson regime [l5| (T ^> d), 
espec ially in nuclei [16j |. but also in several other sys- 
tems [T343 ■ We are not aware of similarly extended and 
precise tests of the RMT approach to chaotic scattering 
in the regime of weakly overlapping resonances (r ~ d). 
Our work is intended to fill that gap. T- invariance vio- 
lation was tested in nuclear spectra (2Cf and for the Eric- 
son regime in compound- nuclear reactions (21rl25T |. Up- 
per bounds on the strength of the T-invariance- violating 
interaction were deduced in both cases. T-invariance vio- 
lation caused by an external magnetic field has also been 
studied in electron transport through quantum dots (26j 
and other devic es | 27l| and in ultrasound transmission in 
rotational flo ws |28| . The RMT approach to T-invariance 
violation [24|, l26l| used in some of these papers is likewise 
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tested very precisely in the present paper. 

The theoretical approach to chaotic scattering is based 
on an expression for the S-matrix originally derived in 
the context of nuclear physics [293 ■ That expression 
contains explicitly the Hamiltonian matrix of the sys- 
tem. Replacing the actual Hamiltonian by a T-invariant 
random-matrix ensemble, one generates an ensemble of 
S-matrices which describes generic features of chaotic 
scattering. Analytical expressions for the S-matrix corre- 
lation functions of that ensemble which apply for all val- 
ues of r J d have been derived [301 ] . These are used in our 
analysis. Replacing the Hamiltonian by an ensemble of 
random matrices with partially broken T invariance [3l| , 
one similarly generates an ensemble of S-matrices that 
describes generic features of chaotic scattering with bro- 
ken T invariance. Some properties of that ensemble have 
been worked out previously 26, 32]. To compare with 
our data we had to extend the theoretical results. This 
work is also reported in the present paper. 

The paper is organized as follows. In Sec. fl] we de- 
scribe the experimental setup and some typical results. 
In Sec. |H] we define the statistical measures in terms of 
S-matrix elements and use these to analyze the data. 
In particular, we define a measure that quantifies the 
strength of T- invariance violation. In Sec. IHII and in the 
Appendix we sketch the derivation of analytic expres- 
sions for the statistical measures. We use the method of 
Ref. [26j]. The theory contains a number of parameters. 
These are fitted to data. We test the theory with the 
help of a goodness-of-fit (GOF) test in Sec. HVDl The 
basic assumption for the applicability of the GOF test is 
that the distribution of the Fourier-transformed S'-matrix 
elements are Gaussian and uncorrelated. In Sec. IIVI we 
demonstrate the validity of that assumption for chaotic 
scattering systems. In the case of T-invariance viola- 
tion, we test the theory further by comparing experi- 
mental values for the elastic enhancement factor and for 
the distribution of the diagonal S'-matrix elements with 
theoretical predictions based on parameter fits to other 
observables. 



I. EXPERIMENT 

For an experimental study of universal fluctuation 
properties of chaotic scattering systems we used flat, 
cylindrical microwave resonators with the emitting and 
receiving antennas acting as single scattering channels. 
As long as the excitation frequency / is chosen below 
/max = Co/ (2 h) , where h is the height of the resonator 
and Co is the speed of light, only transverse magnetic 
TMo modes can be excited, and the electrical field vec- 
tor is perpendicular to the top and the bottom plates of 
the resonator. Then, the associated Helmholtz equation 
for the electric field strength is scalar and mathematically 
identical to the two-dimensional Schrodinger equation of 
a particle elastically reflected by the contour of the mi- 
crowave resonator, i.e., of a quantum billiard [1]. The 



experiments were performed with resonators whose con- 
tour has the shape of a tilted stadium billiard 33] . That 
shape was chosen to avoid bouncing-ball orbits. Each 
microwave resonator was constructed from three metal- 
lic plates. The bottom and the top of the resonator are 
formed by two 5 mm thick high-purity copper plates. 
The center plate had a hole in the shape of a tilted quar- 
ter stadium. The thickness of that plate determined the 
height of the resonator and differed for the two experi- 
ments. The quality factor Q increases with the height 
of the resonator. Thus, to ensure a high Q value the 
center plate for the T-invariant case had a thickness of 
14.6 mm (so that f max = 10.3 GHz). The plate was 
made of brass 12j to technically permit the cutting-out 
of the hole. For the case with broken T invariance a 
copper plate with a thickness of 5.0 mm, which coin- 
cided with the height of the ferrite described below, was 
used (so that / max = 30.0 GHz), see Ref. [HI]. In or- 
der to make sure that only TM modes are excited, the 
excitation frequency / was actually chosen < 25 GHz. 
Screws through the top, middle and bottom plates en- 
sured the good electrical contact needed to achieve high- 
quality values of the resonator. Two thin wires (diameter 
about 0.5 mm) intrude 2.5 mm into the cavity through 
small holes (diameter about 2 mm) drilled into the lid 
of the resonator. They act as dipole antennas to cou- 
ple the rf power into and out of the resonator. A vector 
network analyzer (VNA) provided the rf signal at a vari- 
able frequency / and recorded the signal received at the 
same (other) antenna for reflection (transmission) mea- 
surements. The two signals were compared by the VNA 
in amplitude and phase to determine the complex- valued 
S-matrix elements. These formed the data set for our 
analysis. The cavity is schematically shown in Fig. [T] 
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FIG. 1: The tilted stadium billiard (schematic). The two 
antennas 1, 2 connect the resonator to the VNA. Optionally 
a ferrite is inserted at a fixed location to violate T invariance 
and/or a movable scatterer is used to gather independent data 
sets (see main text). Taken from Ref. [l3 ]. 

For a precise experimental determination of the ele- 
ments of the S-matrix all systematic and statistical er- 
rors must be minimized. The coaxial lines connecting 
the VNA with the cavity are the dominant source for 
systematic errors. They attenuate and reflect the rf sig- 
nal. Both effects were removed by a proper calibration 
of the VNA. Systematic errors are also caused by the 
transmission properties of the two antennas. To account 
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for these, the reflection spectrum of a small cylindrical 
resonator (diameter 5 mm, depth 20 mm) was measured 
using the same antenna geometry as in the actual exper- 
iment. The first resonance is located well above 30 GHz. 
Thus, in the frequency range of interest and for an ideal 
coupling of the antennas to the resonator all rf power 
would be reflected. Any deviation from this expectation 
was attributed to the antennas. The resulting correc- 
tion was applied to the measured spectra in the actual 
experiments. The corrected values of the reflection and 
transmission spectra with the two antennas 1 and 2 pro- 
vided the elements Sn, S12, S21 and S22 of the complex 
2x2 S-matrix as functions of the frequency /. The fre- 
quency step size A/ was > 100 kHz. Typical measured 
reflection and transmission spectra are shown in Fig. [2] 
For the measurement of the S-matrix element Sn(/) an- 
tenna 1 was used as emitting and receiving antenna, for 
that of Si2(/) antenna 2 was used as emitting, antenna 
1 as receiving antenna, etc.. Figure [5] shows that at low 
excitation frequencies the resonances of the billiard are 
isolated, i.e. the mean resonance width T and the correla- 
tion width r, is small compared to the mean level spacing 
d. Since it is a difficult if not impossible task to deter- 
mine the resonance widths in the regime of overlapping 
resonances whereas the correlation width T can be well 
estimated from the data using the Weisskopf formula (36j 
given in Eq. (|28|l below (for more details see Sec. IIIIB[) . 
we refer to the latter in the following. 

As / increases, so does the ratio T/d, and the reso- 
nances begin to overlap. 
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FIG. 2: Reflection spectra (upper panels) and transmission 
spectra (lower panels) of the T-invariant billiard taken at 
three frequency ranges (panel columns). While the left and 
center panels show data at T/d « 0.02 and 0.29, respectively, 
for the two-dimensional regime (where the billiard mimicks a 
quantum billiard), the data at T/d w 0.6 in the right panels 
are obtained in a frequency range where the cavity supports 
three-dimensional field distributions. 

The statistical errors of a single measurement caused 
by thermal fluctuations were reduced by an internal av- 
eraging routine of the VNA. The resulting errors were 
several orders of magnitude smaller than the signal and, 



thus, negligible. The data were analyzed in frequency 
intervals of 1 GHz length yielding M « 10 4 data points 
each. The limited number of statistically independent 
data points in every such frequency interval causes finitc- 
range-of-data (FRD) errors. To increase the number of 
data points and to reduce the FRD errors [H], HH, in 
some of the experiments a small scatterer (an iron disc, 
20 mm diameter) was introduced into the microwave res- 
onator (see Fig. [J) and moved to six different positions. 
We then speak of different realizations of the scattering 
system. 

Experiments with violated Tinvariance were done 
with a magnetized ferrite embedded within the res- 
onator. Such induced Time-Reversal-Invariance Viola- 
tion (TRIV) has been studied in numerous works [l3l . 
|37H42| |. The ferrite has a cylindrical shape (4 mm 
diameter, 5 mm height), a saturation magnetization 
4ttM s = 1859 Oe, and a linewidth AH = 17.5 Oe, 
with 1 Oe=1000/47r A/m. It was provided by courtesy 
of AFT Materials GmbH (Backnang, Germany). Two 
NdFeB magnets (cylindrical shape, 20 mm diameter and 
10 mm height) were placed outside the billiard at the 
position of the ferrite to provide the required magnetic 
fields perpendicular to the top and bottom plates of the 
resonator. The distance between the magnets and the 
ferrite could be adjusted by a screw thread mechanism, 
and field strengths of up to 360 mT could be achieved 
at the position of the ferrite. With this setup TRIV is 
induced via the following mechanism. Because of the 
external magnetic field B the ferrite effectively acquires 
a macroscopic magnetization M that precesses with the 
Larmor frequency ujq around B. This is the origin of 
the ferromagnetic resonance. The rf magnetic field in- 
side the cavity is elliptically polarized and can be de- 
composed into two field components of opposite circular 
polarization with, in general, different amplitudes. Due 
to the Larmor precession of the magnetization M the 
spins of the ferrite couple differently to the two magnetic 
field components. A reversal of time, simulated by an 
interchange of the input and output channels, swaps the 
rotational sense of the two field components and thus, 
due to their different amplitudes, effectively changes the 
coupling of the ferrite to the resonator mode. The in- 
duced TRIV is strongest if the frequency / is close to 
that of the ferromagnetic resonance. The experiments 
with the embedded ferrite demanded a reduction of the 
height of the resonator to 5.0 mm, as the ferrite itself was 
only 5 mm in height. 

In a T-invariant system, the scattering matrix is sym- 
metric, Si 2 = S21. We refer to that property as reci- 
procity. Violation of reciprocity is the hallmark of TRIV. 
In the setup without ferrite the transmission spectrum 
for S2i(/) is indistinguishable from that for S^if) and 
reciprocity holds within the limits given by thermal noise 
(see left panels of Fig. [3j) . Typical transmission spectra 
of the billiard with ferrite and an external magnetic field 
of B = 190 mT are shown in the right panel Fig. [3] The 
two graphs in the upper panels correspond to IS12I 2 and 
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|S2i| 2 , in the lower panel their difference is shown. Fig- 
ure [3] demonstrates that reciprocity is violated. 



II. STATISTICAL MEASURES 

In the present Section we define the statistical mea- 
sures and use them to analyze the data. As pointed out 
in the Introduction, it is our aim to use the data for 
a detailed and accurate test of random-matrix theory. 
Our measures are tailored to this objective. They do not 
address properties of individual resonances but instead 
correlation properties of the fluctuating part of 5-matrix 
elements S a b{f) where a and b take either of the values 
1 and 2. 



A. S-Matrix Correlation Functions 

We decompose the frequency-dependent 5-matrix into 
an average and a fluctuating part, 
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FIG. 4: Three autocorrelation functions for the T-invariant 
billiard determined from the measured values for 5*i2 in the 
frequency intervals 3-4 GHz (T/d w 0.02, circles), 9-10 GHz 
(F/d sa 0.29, triangles) and 14-15 GHz (T/d » 0.6, squares). 
The frequency difference e is plotted in units of the local mean 
level spacing d as obtained from the Weyl formula [4^|. All 
curves are normalized to unity at e = 0. 
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Here and in what follows, the angular brackets (...) de- 
note an average over a suitable frequency interval. In 
order to ensure a more or less constant coupling of the 
electric field modes to the antennas and to the walls of 
the resonator we have always used intervals of 1 GHz 
length. 
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FIG. 3: Spectra of the tilted stadium billiard with a TRIV 
ferrite magnetized through an external field of B = 190 mT 
(left panels) and of B — mT (right panels). The upper 
two panels show the squared modulus of the elements of the 
scattering matrix S12 (full line) and S21 (dashed line) taken 
in the range 15.5-16.0 GHz, where T/d w 0.50. In this range 
the resonator supports, due to its height of 5 mm, only two- 
dimensional modes. Reciprocity is violated for nonvanishing 
magnetic field since S12 7^ 5*2i. To clarify this, we show in 
the lower panels the difference of the squared modulus of the 
elements of the scattering matrix S12 and &i- 



The autocorrelation function of S a b{f) is defined by 

Cab(e) = (S ab (f) S: b (f + e))-\(S ab (f))\ 2 

= (S* b (f)S* b *(f + e)). (2) 

This function quantifies the correlation between S^ b and 
S^b* at two different frequencies / and f + s. Figure 0] 
shows three examples of autocorrelation functions all ob- 
tained from data for the billiard without ferrite. The rate 
of decrease of the functions with increasing e depends on 
the ratio T/d. None of the functions has the Lorentzian 
shape predicted by Ericson [n| for the regime of strongly 
overlapping resonances T ^> d. We show later that the 
rate of decrease agrees with random-matrix predictions 
for the relevant values of T/d. 

To quantify TRIV we measure the violation of reci- 
procity by the cross-correlation function of 5i2(/) and 
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((Sf 2 (f)S^*(f + e)) 

V(\sf 2 (f)\ 2 )(\si(fW) 



(3) 



For a T-invariant system, reciprocity holds, and 
Ccross(O) = 1. In case of complete TRIV we expect that 
5i2 and 52 1 are completely uncorrelated, C cr0 ss(e) = 
for all values of e. (This expectation is borne out in 
Sec. IIIIl see also Ref. (42| for a treatment of the two-level 
case). In summary we have 



Ccross (0) — 



1 for Tinvariance , 
for complete TRIV 



(4) 



As explained in Sec. HI we have used six realizations to 
increase the statistical significance of the data. For each 
realization the cross-correlation coefficient C cross (0) was 
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FIG. 5: Cross-correlation coefficient C cr oss(0) as a measure of 
TRIV. In each frequency interval of 1 GHz length C crO ss(0) 
was evaluated as an average over six realizations. The points 
denote the mean values and the error bars denote the standard 
deviations. The zero on the ordinate is suppressed. Based on 
Ref. [H|]. 



in Fig. [6] show a fit of the random- matrix expression de- 
fined in Sec. lIIII to the data (the fit procedure is described 
in Sec. IIV Ap and correspond to the local-in-time mean 
values of the Fourier coefficients. The data are seen to 
scatter about their time-dependent mean. In Sec. IIV Bl 
it is shown that the data points divided by their local 
mean value at different times are indeed uncorrelated and 
that the distribution of the rescaled Fourier coefficients 
of the autocorrelation function is exponential. The decay 
of the average function (solid line in Fig. ^ is faster for 
T/d = 0.29 (lower panel, frequency interval 9-10 GHz) 
than for T/d = 0.02 (upper panel, frequency interval 3- 

4 GHz). This is due to stronger absorptive losses. In 
both cases the decay is non-exponential (and the auto- 
correlation function is, therefore, not Lorentzian). At 
T/d = 0.29 and for times larger than about 1000 ns the 
decay is dominated by noise. Nevertheless, a decay over 

5 orders of magnitude is experimentally well established. 



computed and the average over all realizations was taken. 
For an external magnetic field of 190 mT the resulting 
averaged cross-correlation coefficient is shown in Fig. [5] 
This coefficient deviates noticeably from unity around 6, 
16 and 24 GHz, indicating TRIV. The first dip can be at- 
tributed to the ferromagnetic resonance which in our case 
is located at 6.6 GHz. We assume that the other two dips 
arise from an enhancement of the influence of the ferro- 
magnetic resonance by standing rf magnetic fields inside 
the ferrite. Yet, the smallest values (C cross (0) ~ 0.4) ob- 
tained are well above zero. Hence, at a field strength 
of 190 mT the ferrite induces only a partial violation of 
T invariance. This is also found at the other investigated 
field strengths up to 340 mT. In Sec. IIIII we show that 
the strength of TRIV can be deduced from C cross (0). 



B. Fourier Transformation 

The measured scattering matrix elements are corre- 
lated for neighboring frequencies /. The correlations re- 
sult in a non-zero value of the autocorrelation function 
defined in Eq. ([2]) and depicted in Fig. [4] We show in 
Sec. IIV Bl that after a Fourier transformation the cor- 
relations between data points at different times can be 
removed. This facilitates a statistically sound analy- 
sis and is our motivation for using that transformation. 
Since S a b{f) is measured at a discrete set of frequencies, 
the Fourier coefficients S a b(k) are likewise obtained at 
discrete time points tk — k/A, and the same is true 
of the autocorrelation function C a b(s) and its Fourier 
transform C a b{k). Here, A = 1 GHz is the length of 
the frequency interval and k = 0,1, . . . , M — 1, see Sec- 
tion llV A1 We simplify the notation by using as argument 
of the Fourier transforms the integer k. According to the 
Wiener-Khinchin theorem we have C a b(k) = \S a b(k)\ 2 . 
Figure |5] shows two examples of C\2{k) at different val- 
ues of T/d for the 7~-invariant system. The solid lines 



C. Elastic Enhancement Factor 

In chaotic scattering, elastic processes are known to 
be systematically enhanced over inelastic ones. The ef- 
fect was first found in nuclear physics [44j-|46| but plays 
a role also in mesoscopic physics 27]. The enhancement 
depends on the degree of T violation. The elastic en- 
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FIG. 6: Fourier coefficients of the autocorrelation function of 
512 (in semi-logarithmic scale). Data points are from the bil- 
liard without ferrite in the frequency ranges 3-4 GHz (upper 
panel) and 9-10 GHz (lower panel). For clarity only every 5 th 
data point is shown. The solid lines are best fits to the data. 
In the data shown in the lower panel the decay is dominated 
by noise for times larger than about 800 ns. These data are 
not taken into account in the fitting procedure. 
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hancement factor is defined as 



\\ 2 ) 



W = sj{\^?) (\S^)/(\Sf 
= VCii(0)C 2 2(0)/C , i2(0), 



(5) 



where the second equality results from Eq. In the 
limits of isolated resonances with many weakly coupled 
open channels and of strongly overlapping resonances the 
values for W are 



W 



1 + 2//3 forr/rf<l 
2/f3 forr/d»l. 



(6) 



Here, /3 = 1 for T-invariant systems and (3 — 2 for com- 
plete TRIV. The elastic enhancement factor VV was de- 
termined in two ways: (i) Using the first of Eqs. ([3]) we 
calculated the averages over frequency directly from the 
experimental values for S a b(f)- This amounts to deter- 
mine W from a single experimental value for each of the 
autocorrelation functions Cn(0), Ci2(0), 6*22(0). (ii) In 
the second of Eqs. ([5]) we used the values of the autocor- 
relation functions obtained by a best fit of the analytical 
expression given in Eq. (1261) below to the experimental 
one. These are the solid lines in Fig. |S1 The method 
of fit (described in Sec. IIV A|) uses the entire data set 
and is, therefore, expected to give more reliable values 
for W. This is indeed borne out by the results shown in 
Figs. [7J and |U For the T-invariant case shown in Fig. [7J 
the elastic enhancement factor decreases from W f=s 3 at 
low frequencies (r <C d) to W ~ 2 at high frequencies 
(f « d), in qualitative agreement with Eq. ©. 
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FIG. 7: Elastic enhancement factors W for the T-invariant 
billiard. The open circles are obtained with method (i), the 
filled circles with method (ii) described in the text. The error 
bars indicate uncertainties due to the finite range of data [34l . 
[35I ]. Above 10 GHz the analogy to a quantum billiard breaks 
down and Eq. ()26|) needed for the analytic evaluation of W is 
no longer applicable. The dashed horizontal lines indicate the 
limits of W for T-invariant systems: Upper line for F <C d, 
lower line for r ^> d. 



from 3 to 2 with increasing T/d. However, at frequencies 
of about 6, 16 and 24 GHz dips are observed. Around 
16 and 24 GHz the values of W drop below 2. This is 
not possible for a T-invariant system. These features 
are similar to those of the cross-correlation coefficient in 
Fig. (3) Both measures indicate a substantial violation of 
T-invariance at about 6, 16 and 24 GHz. 
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FIG. 8: Elastic enhancement factors W for the billiard with 
partial TRIV. The open circles are obtained with method (i), 
the filled circles with method (ii) as described in the text. 
The error bars show the root-mean-square values for the 6 
realizations. The dashed horizontal lines mark the limits of 
W for T-invariant systems as in Fig. [7J Taken from Ref. 



III. THEORY 

As stated in the Introduction, it is the aim of the 
experiments reported and analyzed in this paper to 
test random-matrix theory as applied to chaotic scat- 
tering systems. According to the Bohigas-Giannoni- 
Schmit conjecture Q , the spectral fluctuation properties 
of chaotic T-invariant quantum systems coincide with 
those of the Gaussian orthogonal ensemble [GOE], those 
of quantum systems with complete TRIV with those of 
the Gaussian unitary ensemble [GUE] of random ma- 
trices. Systems with partial violation of Tinvariancc 
are accordingly described by a crossover from orthogo- 
nal to unitary symmetry. For such systems, analytical 
expressions for central statistical measures of scattering 
processes (i.e., the autocorrelation function @ and the 
cross-correlation coefficient ([3])) have not been worked 
out before. We fill that gap in Sec. IIII Al 



A. Crossover from Orthogonal to Unitary 
Symmetry: the Autocorrelation function and the 
Cross-Correlation Coefficient 



Results for the billiard with violated Tinvariance are 
shown in Fig. [8] Although W was obtained from a data 
set of 6 realizations, the values obtained with method 
(i) still show large uncertainties while method (ii) yields 
reliable results. Again W displays an overall decrease 



For chaotic scattering processes, the GOE — > GUE 
crossover was extensively investigated in Ref. [26| ■ There 
the destruction of weak localization by an external mag- 
netic field in the transmission of electrons through a few- 
channel disordered microstructure was determined. The 
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connection between the conductance g and the S-matrix 
is given by the Landauer formula, 



A/2 A 

E E 

a=l 6=A/2+l 



l^a| 2 } 



(7) 



where A counts the total number of open channels and 
an equal number of incoming and outgoing channels ^ is 
assumed. Efetov's supersymmetry method [49j was used 
to calculate ensemble averages of squares of S-matrix el- 
ements S a b for a ^ b. The S-matrix embodies the GOE 
— > GUE crossover in the manner described below. Here 
we describe the extension of that approach to the calcu- 
lation of the autocorrelation function @ and the cross- 
correlation coefficient ([3]). These observables were not 
considered in Ref. [26|. 

We write the unitary scattering matrix in the general 
form [H 



N 



Sab(f) = S ab - 2wi J2 W W [ D ^} uv W b> 



(8) 



where the inverse propagator D is given by 



c=l 



Wc, 



(9) 



The matrix elements W a ^ and W oll describe the coupling 
of antennas a and b with the resonator mode pt [42( . The 
sum over c in Eq. ([9]) extends over the two antennas but 
includes also a number (A — 2) of fictitious channels. The 
latter describe Ohmic absorption in the cavity |l9[ . The 
matrix H^ v is the Hamiltonian of the closed billiard. It 
has dimension N and the limit N — > oo is eventually 
taken. 

The coupling matrix elements W cll are chosen real, 
W c ,j, = W* for all c,fi, and a violation of Tinvariance 
by an external magnetic field is taken into account only 
in H^ u . We also assume that the W cli are independent 
of frequency /. That assumption holds within frequency 
intervals of 1 GHz width. The experiment was designed 
such that direct power transmission between the anten- 
nas is excluded so that the average S-matrix is diagonal. 
We have verified that fact experimentally. A diagonal 
average S'-matrix is implied by the relation 



JV 

E 



W a ^W b ^ = Nv 2 a 5 ab . 



(10) 



The parameter v\ measures the average strength of the 
coupling of the resonances to channel a. The Hamiltonian 
H^ v is a member of a random-matrix ensemble describ- 
ing partial violation of Tinvariance. In random-matrix 
theory, the GOE — > GUE crossover is written as [3l| 



11 ^v — 11 fit, i l J ' u 



(11) 



The real and symmetric matrix is a member of the 
GOE, and the elements of the real and antisymmetric 
matrix are uncorrelated Gaussian-distributed ran- 

dom variables. Thus, 

(H$) = «»} = , 

= ^(V<W+V<V) - 

(H$>H$,) = ^0V<W- W) ■ (12) 

Here A has the dimension energy and for the GOE de- 
notes half the radius of Wigner's semicircle. The param- 
eter £ measures the strength of T-invariancc violation. 
For ■K^/yN = 1 the matrix H is a member of the GUE. 
However, on the local level (energy intervals measured 
in units of the mean level spacing d of the GOE) the 
transition from GOE to GUE already takes place when 
the typical matrix element of the TRIV term becomes 
comparable to d — ttX/N, i.e. when 



A 



N VN 



N VN 



(13) 



or when £ ~ 1. The S'-matrix (JHJ is symmetric only for 
£ = 0, and reciprocity does not hold for £ ^ 0. 

Starting with Eq. (QJ we have used angular brackets 
to denote the running average over (parts of) the exper- 
imental spectra. The averages were actually taken over 
1 GHz frequency intervals. Now we consider an ensemble 
of S-matrix elements of the form of Eq. ([8j obtained by 
inserting many realizations of the random Hamiltonian 
H^ u in Eq. ([TT|) . We denote averages over that ensemble 
also by angular brackets. This is legitimate because er- 
godicity guarantees the equality of ensemble average and 
of the running average over a single realization of the 
ensemble, see Sec. II C3 of [47j . 

The autocorrelation function for the S-matrix defined 
in Eq. (|5J) is known for the case of Tinvariance (£ = 
0) [H, for the case of complete TRIV (£ = vW/tt) [H 
and had to be calculated for the case of partial TRIV 
(0 < £ < yN /it). For the first case it reads 



1 



OO pOO 



o Jo Jo 



cxp (~*"^-(a*i + A*2 + 2/x) 

TT l-Tc/J, 

11 v/a+^MiXi+T^) 



x J ab (fi, nx, fi 2 ) ■ 



(14) 



The upper index of C a b indicates that an average over 
the GOE was taken. The integration measure is given by 



M(l ~ - Mgj 
(M + M1HM + M2) 2 
1 

\/(Mi + M1W1 + M2) ' 



(15) 
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and we have 
Jab(^, ^1,^2) 



<U 1 (S aa ) I 2 T 2 




( Ml + 


(12 , 2/x \ 


\l + T atH 1 


+ T a H2 I — T a il) 


(1 + 8 ab ) T a T b 




.(i + T B ^i)(i + r 6 /*i) 


^2(1 + ^2 





(i + T a/ i 2 )(i + r bM2 ) 

2/i(l - A*) 



(l-T a /i)(l-T 6 /t) 



(16) 



Here and in Eqs. (1191) and (|24l) below the input param- 
eters are the mean level spacing d and the transmission 
coefficients T c in all channels c, defined as 



T 



1 



(ft 



(17) 



We observe that < T c < 1. For all three cases 
the matrix elements W c/1 occur in the final expression 
for the correlation functions only via the transmission 
coefficients. Using the analytical result for (S aa ) = 
(1 - Tr 2 v 2 c /d)/(l + n 2 v 2 /d) and Eq. (HUJ) one finds 
that 



T 



4ir 2 v 2 /d 

(l+W 2 V 2 /d) 



(18) 



The choice of the parameters T c is described in Sec- 
tion IIV Al The threefold integrals in Eq. (fT4|) and in 
Eq. below are numerically computed most conve- 
niently in terms of the integration variables introduced 
in Ref. go|. 

For the second case, the ^-matrix autocorrelation func- 
tion was worked out in Ref. [HI, 



°a6 



(0 



00 />1 

d/ii / d/iexp 
Jo 



.2ire. . 
j-(Mi +N 



ni 



1 - T c fi 



T a 



T b 



{l + Taflx^l-Tafl) (1 + TbtnXl-Tbfl) 



1 



Mi + A* 

/ii/i(T Q + rf,-r a Tf,)} 



{mi - /* + 1 



(19) 



The upper index of C a b now indicates the average over 
the GUE. 

In our experiments we deal with partial TRIV, i.e. 
with the third case and the iS-matrix autocorrelation 
function had to be calculated for all values of the param- 
eter £ introduced in Eq. (TlTj) . To this end we generalized 
the work of Ref. 26|. We present here only the result 
and defer details to the Appendix. The autocorrelation 
function is given in terms of a threefold integral over inte- 
gration variables Aq,Ai, and A2, see Eq. (2) of Ref. [5lj |. 



However, the integrals are evaluated numerically more 
conveniently in terms of the integration variables given 
in Sec. 5 of Ref. [5(| ■ For the transformation to these 
one needs to distinguish in the integrations over Ai and 
A2 the case where Ai > A2 and the case where Ai < A 2 . 
For instance, for the case Ai > A2 the transformation to 
integration variables fi, fix, fi2 is given by 

A = l-2/i, 

Ai = V(l + A*i)(l + A*2)+MiM2+W, 



where 



A 2 = v / (l + Ati)(l + At2)+/ii/i 2 -W, (20) 



U = 2VMi(l + AtiWl + /t 2 ). 



Then both, the S'-matrix autocorrelation function and 
the cross-correlation coefficient are obtained as special 
cases of a function F£ b (e). With the notations 



and 



t = 7T 4 , 



K = 4(/x + /ii)(/t + /i 2 ) , 

t = 4/x(i - /t), g = a 2 - 1, n = a 2 - 1, 

e± = 1 ± exp(-2tF), 
7 _ (2-T a )A 2 +T a Ai 

A a 
Ba 
Ca 
C 2 



(21) 



(22) 



4(l + T a Ai)(l + T a A 2 )' 

(2-r a )Ai+r Q A 2 
4(l + T a Ai)(l + T a A 2 )' 

1 1 r = m(i-m) 

21-T /i' ' 
W / 1 



(1 - T at i) (1 - Tfc/t) ' 
1 



4 Vl+Ta/i2 l + Tfc/t! 
1 1 

+ - 



1 + T a /ti 1 + T 6 /i 2 / ' 
the function F° b (e) reads 



(23) 



d/ii / d/t 2 / dfi?^ 1 ' 1 ' 1 ' 1 '" 



11 Jo Jo 

iire 



T 



x exp 



'-(fll + /t 2 + 2/t) 
1 - T c /i 



VCl + Tc/iOa + rc/ia) 
x [exp(-2tW) • {j^0i,A*i,A*a) • [Fe+ 
+ (A 2 - A 2 )e_ + 4t^(A 2 e_ + F{e+ - 1))] 
+ a ■ 2(1 - <U)T a T b # o6 } + (Ai o A 2 )] . (24) 
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Here, 



K ab = 



e_ 2iUC 1 T (25) 
2jr{(i C b + A b C a )Q\2 + {B a C b + B b C a )H\i} 
3d J 7 - C 2 {\\ - Af) + C 2 m(4X 2 2 - IT) 
(e+ - £ —) [3d(A^ - Af) + tf«7i(4A2 - 2 J 7 ) 
2J-{(I a C 6 +i 6 (7 a )£A 2 - (B a C b + B b C a )UXi] 
(2tK - 1)C 2 T ■ 



The integration measure J^(fJ,, fj,%, fj, 2 ) arL d the func- 
tion J a b(n, Mi, M2) ar e given explicitly in Eqs. (|T5|) and 
Eq. (HU). Setting er = + (cr = -) in Eq. (jUJ yields 
the autocorrelation function (the cross-correlation coeffi- 
cient), 



C ab (e) = F+{e) 
,(e = 0) = F" 6 (e = 0) 



(26) 
(27) 



We observe that for £ = 0, i.e. t = 0, the function K ab 
defined in Eq. (|23|) vanishes and F+ b (e) in Eq. (j2"I)) turns 
into the autocorrelation function of the GOE given in 
Eq. (|l"4|). We checked our analytical results by compari- 
son with RMT simulations. In Fig. [9] we show the cross- 
correlation coefficient versus £ as obtained analytically 
and by RMT simulation for a typical set of transmission 
coefficients. We also indicate how the analytical result is 
used to determine the value of £ from a measured value 
of the cross-correlation coefficient. To test the validity of 
Eq. (f2~4")) we compare in Figs. [TU] and [IT] analytic results 
for the autocorrelation functions with numerical simula- 
tions, both in the frequency and in the time domains. 
The parameter r a b s measures absorption and is defined 
in Sec. IIII Bl below. In all cases, the agreement is very 
good. 




FIG. 9: Dependence of the cross-correlation coefficient 
Ccross (0) on the parameter £ as predicted by a random-matrix 
model for partial violation of Tinvariance. The analytic re- 
sult (line) is compared with an RMT simulation (dots) for the 
same set of transmission coefficients. Also shown is how an 
experimental value of C cr oss(0) = 0.49(3), c.f. Fig. [5] trans- 
lates into £ = 0.29(2). Based on Ref. Qj. 



The theoretical expressions given in Eqs. (fTTj) . (|19|) and 
([24| are obtained as averages over the ensemble of Hamil- 
tonian matrices defined in Eq. (|III) in the limit N — > 00 




(MHz) 



FIG. 10: Comparison of the analytic result for the autocor- 
relation function C11 versus e for transmission coefficients 
Ti = 0.407, T 2 = 0.346, r abs = 2.41 and TRIV parameter 
^ = 0.293 (solid line) with RMT simulations (dots). We show 
only the result for Cn as that for C12 is barely distinguish- 
able. The curve is normalized such that it equals unity for 
£ = 0. 
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FIG. 11: Comparison of the analytic results for the Fourier 
transform of the autocorrelation function versus time t for 
transmission coefficients Ti — 0.407, T2 = 0.346, r a b s = 2.41 
and TRIV parameter £ = 0.293 to RMT simulations. We 
show the results for Cat for a = b = 1 (dashed line and filled 
points) and a = 1, b — 2 (solid line and crosses, respectively). 



and directly yield the autocorrelation function. In con- 
trast, Fig. [4] shows autocorrelation functions obtained by 
averaging the data in a frequency interval of 1 GHz width. 
For conceptual clarity we distinguish both cases by refer- 
ring to the theoretical and to the experimental autocor- 
relation functions, respectively. 



B. Parameters 



The parameters in Eqs. yjj, UB and C31> are tne 
average level spacing d, the transmission coefficients T c 
for all channels c, and the parameter £ for TRIV. We 
have calculated d from the Weyl formula 43j . A starting 
value for the parameter £ was determined from the ex- 
perimental cross-correlation coefficients shown in Fig. [5] 
as described in the caption of Fig. |H1 Here we use that 
the cross-correlation coefficient depends only weakly on 
the transmission coefficients in the frequency range 1- 
25 GHz. Results are shown in Fig. [T2] The largest value 
of £ is £ ~ 0.3. In determining T c by fitting the theoret- 
ical expressions for the autocorrelation function to the 
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FIG. 12: Values of the TRIV parameter £ for the billiard 
with the ferrite magnetized with B = 190 mT. The error bars 
indicate the variability of the results within the 6 realizations. 
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data, we also use £ as fit parameter, with starting value 
as just described. For the channels c = 1 and c = 2, 
i.e. for the antennas, we also have determined starting 
values from the definition Eq. (|17p and from the mea- 
sured values of Sii(/) and 822(f)- The fits discussed in 
Sec. QV] yielded T\, T 2 within 5 % of these starting values. 
The remaining transmission coefficients describe Ohmic 
absorption in the walls of the resonator and the ferrite. 
If T c <C 1 for all absorbing channels, the products over 
these channels appearing in Eqs. ([T3|l . f(T9]) and (|24|) sim- 
plify so that each of the three theoretical autocorrelation 
functions depends only on the sum r a b s of the transmis- 
sion coefficients for the absorbing channels. Accordingly, 
in addition to T\ , T 2 and £ the parameter r a b s was used as 
fitting parameter. To estimate the correlation width Y of 
the resonances, we have used the Weisskopf formula [36[ 



2 7T - = ^ T c = Ti + T 2 + T a b s 



(28) 



and the fitted values for T\, T 2 , and T a b s . Using nu- 
merical simulations and with help of the analytic result 
Eq. (fl4|) for the autocorrelation function we checked that 
this formula indeed yields a very good estimate for the 
correlation length even in the regime of weakly overlap- 
ping resonances and for a few open channels. The re- 
sults of our fits (each done in a frequency interval of 1 
GHz length) are displayed in Figs. H3] and Q31 We show 
Y/d as obtained from Eq. (|28|) (top panel), r a b s (middle 
panel) and Ti and T 2 (bottom panel) versus frequency 
for the case without and with TRIV, respectively. The 
transmission coefficients and r a b s generally increase with 
frequency. We note that without the ferrite Y/d never 
exceeds the value 0.3 (regime of weakly overlapping res- 
onances) since the excitation frequency must be chosen 
below / max = 10.3 GHz, whereas in the system with fer- 
rite we have / max = 30 GHz and thus Y/d attains values 
as large as 1.2. 



C. Distribution of ^-Matrix Elements 

The distribution of S'-matrix elements is not known 
completely in analytical form, neither for the T-invariant 



FIG. 13: Y/d (top panel), r a b s (middle panel) and the 
transmission coefficients T\ (filled circles) and T2 (open cir- 
cles) (bottom panel) versus frequency for the billiard with 
Tinvariance. The errors are typically of the size of the sym- 
bols. 
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FIG. 14: The same as in Fig. [13] but for the billiard with the 
ferrite, obtained as averages over six realizations. The scatter 
of the values for different realizations about the mean value 
is typically of the order of the symbol size. 



system nor for TRIV. The most complete information is 
available in the regime of strongly overlapping resonances 
(r/d > 1). In a basis in channel space where (S a b) is 
diagonal, the inelastic S'-matrix elements S a b = S^ b with 
a ^ b have a bivariate Gaussian distribution (52j. Thus 
the phase of S a b is uniformly distributed in the interval 
{— 7r,7r}. The distribution function P of the modulus 
r = \S a b\ depends only on the ratio z = r/((|S' a b| 2 )) 1 / 2 
and is given by 



P(z) 



-z exp 



(29) 



The diagonal elements S aa have a bivariate Gaussian dis- 
tribution only for r ^ d and only if K^a)! <C 1. Other- 
wise, unitarity constraints cause the distribution to differ 
from the Gaussian form [o3l - [55l |. 

In Ref. [13] an analytic expression for the distribution 
of the elastic elements S aa of the scattering matrix of a 
generic chaotic system without or with partially violated 
Tinvariance was derived. It applies for cases with many 
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open channels. With the notation 

■a + 1 



Saa — y/faG a j %a 



V 9a ~ T a X ' 



the distribution P{x a ,9 a ) of S aa is given by 
„ df(») 



(30) 



•(31) 



With the help of the transformation Eq. ([20j) . the defini- 
tions (|2"3")1 and with u> = (j/ — 1)/2 the analytic expression 
for the function f a (y) given in Ref. [32[ may be cast into 
the form 



d^ 2 / d^i / d/i- 



J" 



x exp 



(-^(M!+M 2 + 2 M )) 



w + n 



V( w - - w) n ^(i + r c ^i)(i + r c /i 2 ) 

exp (-2W) • [^6+ + (A2 - A?)e_ + 4tft(A§e_ + T(e+ - 1))] + (Ax <-> A 2 ) 
For 7"-invariant systems the threefold integral can be simplified, and P(x a , 9 a ) takes the form 

P{x a ,9 a ) = J- -T- (1+2/) [ T ab S {K 1 {w)J 2 { W ) + K 2 (w)J 1 {w)) 

4-7T ay L 
A 

+ ^2%(Lt(w)H2(w)+LZ(w)HZ(w)) 



l-T c ji 



c=l 



Here, t a c = 1 for c = a and i° = T c otherwise, A is the number of open channels and 

coo e -r aba y/2 



dy- 



Vv\y~ w \ d=i V 1 + t dy 



dy 



1 f ' - 



Vy\y^\L\v T +tfy :i + ta cy 



dye' 



ab= y/2_ 



vy\y- 



Ud=i y/l + %V [ v 



tII(i-«)-7 



dye" 



ab= y/2_ 



d=£b 

\Zy\y-w\ 



nti xA+^ 



r II( 1 -^)-- 

1 . . w 



d=£b,c 



(32) 



(33) 



(34) 



The corresponding functions with index 2 are given by 
the same expression except that the integration limits 
w, 00 have to be replaced by 0, w. 

These analytic results were previously tested experi- 



mentally for the case of a single open channel plus ab- 
sorption in Refs. [l~8l. l56l - [59| . In Fig. [15] we compare for 
several frequency intervals the experimental distributions 
of the elastic 5-matrix element Sn with the theoretical 
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predictions for the case of two open channels with ab- 
sorption. The data were taken with the billiard used for 
the experiments with partial TRIV but without the fer- 
rite. This was done because that billiard has a smaller 
height so that the range where only a single vertical mode 
is excited, extends up to 30 GHz. Higher values of the 
frequency result in larger absorption and in larger values 
of T/d. It is here that the theoretical result (j3"3")) is ex- 
pected to apply. The value of T/d was determined from 
the Weisskopf formula given in Eq. (|28p and from the 
values of T±, T 2 , r a b s obtained from a fit of the Fourier 
transform of the S-matrix autocorrelation function as de- 
scribed in Sec. IIVI The very good agreement corrobo- 
rates the precision of the fitting procedure and of the 
GOF test discussed below. Note that the distributions 
are very far from a bivariate Gaussian distribution. 

An analytic expression for the distribution of the off- 
diagonal elements of the S-matrix exists in the Ericson 
regime r>d, however not in the range of T/d achieved 
in the experiments. In Fig. [16] we, therefore, compare 
experimental distributions to RMT simulations. We note 
again the good agreement. In the frequency range 23- 
24 GHz (where T/d 1.01) the distribution of |S 12 | is 
well described by Eq. (f2U)) , and the distribution of the 
phases is nearly uniform. Thus, in this frequency range 
the distribution of the non-diagonal S'-matrix elements is 
already close to that expected in the Ericson regime while 
that of the diagonal elements is still far from Gaussian. 
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FIG. 16: Distribution of Si2-matrix elements according to 
modulus (upper panels) and phase (lower panels) for data 
points as described in the caption of Fig. 1151 An RMT simu- 
lation (solid lines) shows acceptable agreement with the data. 



IV. DATA FITS AND DISTRIBUTION OF 
FOURIER COEFFICIENTS 

In the present Section we test predictions of random- 
matrix theory with the experimental data. We proceed 
as follows. Using the results of Sec. IIIB1 we fit the 
parameters of the Fourier transforms of the theoretical 
expressions for the autocorrelation functions (Eqs. (frjj 
and (|2"4"|) ) to the data. We show that within the accu- 
racy of the data and after rescaling, the distribution of 
the Fourier-transformed S'-matrix elements is Gaussian. 
This property is used to develop a goodness-of-fit (GOF) 
test that quantitatively tests the quality of RMT predic- 
tions. As a second test of RMT we compare predicted 
values of the elastic enhancement factors with the data. 

In Sec. IIII CI we have shown that the real and imag- 
inary parts of the S-matrix elements in general do not 
have a Gaussian distribution. How can this fact be rec- 
onciled with the statement just made that the Fourier- 
transformed S-matrix elements do have such a distribu- 
tion? The Fourier transform is a linear transformation, 
after all. We answer that question as we proceed. 



FIG. 15: Distribution of Sii-matrix elements according to 
modulus (upper panels) and phase (lower panels). The his- 
tograms give the probability distribution functions in the fre- 
quency ranges (from left to right) 9-10 GHz, 17-18 GHz and 
23-24 GHz. The data were measured with the billiard used 
for TRIV but without ferrite and for a total of 8 realizations, 
i.e. each graph is constructed from 80 000 data points. All 
of these were used in the histograms, and correlations be- 
tween S'-matrix elements at neighboring energies were thus 
neglected. Nevertheless, the analytical result Eq. (|3ip (solid 
lines) agrees with the data. 



A. Fits 

We focus attention on the fluctuating part S (/) of 
the S-matrix elements (see Eq. ([T])) and omit the indices 
a and b for brevity. By definition we have (S fl (f)) = 0. 
Data are taken at frequency increments A > 100 kHz. 
The mean level spacing d, the transmission coefficients T\ 
and T2, and the absorption coefficient r a b s are typically 
constant in frequency intervals of 1 GHz width. In every 
such interval we have M ~ 10 4 measured values of S fl (/) 
for all combinations of channel indices a, b. We write 
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fj = fa + jA where /o is the frequency at the lower 
end of the interval and j = 0, 1, 2, . . . , (M — 1). We use 
discrete Fourier transformation and define 



M-l 



, M - 1 . 



(35) 



so that 



M-l 



M 2^ 6 



3=0, 



,M-l. (36) 



k=Q 



We fit the parameters in the theoretical expressions (|14p 
and (IM1) to the distribution of the squares Xk = |S\| 2 of 
these Fourier coefficients. The Wiener-Khinchin theorem 
states that the latter are equal to the Fourier coefficients 
C(k) of the experimental autocorrelation functions C(e). 
We accordingly calculate the Fourier transforms C(k) of 
the theoretical autocorrelation functions in Eqs. (|14p and 
([24| at the same discrete values of k as occur in the dis- 
crete Fourier transformation in Eq. (|35p . The parameters 
are the transmission coefficients T±, r a b s and £. For 
the transmission coefficients 7\ and T 2 we used Eq. (fl~7|) 
with experimental values for (S aa ) and a = 1, 2 as start- 
ing points but allowed the values of Ti and T 2 to vary. 
The best-fit values differed by no more than 5 % from the 
starting values. For £ a starting value was obtained from 
a comparison of the experimental and the analytic cross- 
correlation coefficients as outlined in Sec. IIIIBI For the 
fit parameter r a b s no starting values could be computed 
from the measured data. The fit yields the solid lines 
shown in Fig. [5] and defines (xfc). We observe that (xt) 
decreases by several orders of magnitude over the avail- 
able range of k values. We also observe that for k 7^ k' 
the Fourier coefficients S& and Sk> are uncorrelated. This 
follows from the fact that the autocorrelation functions 
depend only on the difference e of the two frequency ar- 
guments and is shown below in Eq. (|39[) . 

As stated in the Introduction, the use of generic expres- 
sions derived from RMT is justified only for energy spac- 
ings bounded from above by the period of the shortest 
periodic orbit in the classical microwave billiard. There- 
fore, Fourier coefficients are generic only for times larger 
than the repetition time of the shortest periodic orbit 
We have estimated that time and found it to be smaller 
than the time-index k — 5 of Sk in Eq. (|35|) . We note, 
however, that the average S-matrix elements (which cor- 
respond to k = 0) are not generic. This is mirrored by 
the fact that these are used as input parameters in our 
analysis. 



B. Gaussian Distribution 

We ask: How are the |S\| 2 distributed about their 
mean values determined by the fits? In order to study the 
distribution of the S\ and of the coefficients Xk = \Sk\ 2 
with good statistics, we must sample all data points 



shown in Fig. [6j To this end we remove the strong and 
systematic k dependence by rescaling: We divide Sk by 
\J (xk) and Xk by (xk) and find that the renormalized 
S'-matrix elements Sk/ \f {xk) have a bivariate Gaussian 
distribution both for the elastic (case shown in Fig. [P7|l 
and the inelastic one (shown in Fig. [T5]) . The left-hand 
side of Fig. Qj|] shows that after rescaling of the Xk the 
logarithms of the rescaled coefficients Zk — Xk/ (xk) scat- 
ter about zero. Moreover we find that the distribution of 
the Zk is stationary in k. By this we mean that the distri- 
bution of the Zk determined from sampling their values 
within some interval of length Sk <C M does not depend 
on the choice or length of that interval. The statistical 
accuracy of that statement is obviously limited by the 
fact that the number of data points contained in the in- 
terval decreases with decreasing length Sk. Stationarity 
allows us to study the joint distribution function of all Zk 
obtained from S'-matrix data that lie within a frequency 
interval of length 1 GHz. That step improves the statis- 
tical accuracy of the result. 
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FIG. 17: Distribution P of the rescaled Fourier coefficients 
5/v/ (xk) for the elastic case, {a, b} = {1, 1}. Upper panels: 
The data (histograms) for the distribution agree well with the 
solid lines given by Eq. (|29p . Lower panels: The phases are 
uniformly distributed in the interval {— 7r,+7r}. Data source 
is as in Fig. [TS] with Fourier coefficients taken from the first 
200 ns, i.e. 1 240 data points contribute to each histogram. 



The right-hand side of Fig. [19] shows that the coeffi- 
cients Zk have an exponential distribution, as expected 
for the absolute squares of variables with a bivariate 
Gaussian distribution. To test this statement quantita- 
tively we observe that for an exponential distribution the 
ratio (zl)/(zk) 2 should have the value two. For our finite 
data set we define 
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FIG. 18: Same as in Fig.[T7]but for the inelastic case {a, b} = 
{1,2}. Data source is as in Fig. \Tft\ with Fourier coefficients 
taken from the first 200 ns. 
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FIG. 19: Distribution of the Fourier coefficients Xk = S12I 2 
in the interval 16-17 GHz under TRIV with B = 190 mT. 
Panel (a) on the left-hand side displays on a logarithmic scale 
(base 10) the x^ as obtained from the data. Panel (b) shows 
analogously the rescaled quantities — Xk/{xk)- Panel (c) 
on the right-hand side shows the distribution of M = 1200 
rescaled coefficients on a logarithmic scale for six realizations. 
The dashed line is the exponential expected for a Gaussian- 
distributed 5*. 



and obtain for the variance of M 2 /Mi 




(38) 



Evaluation of that ratio for the data set in the range 
16-17 GHz shown in Fig. [19] with 6 realizations gives 
M 2 jM\ = 1.97, which is within the defined error limits 
(2 ± 0.09) for the M = 6 • 200 = 1200 contributing data 
points. A systematic analysis of our data ensemble for all 
24 frequency intervals between 1 and 25 GHz yields 13 
accepted and 11 rejected ratios, i.e. 54 % of all frequency 
intervals are within the 1-a range defined by Eq. (|3"8")) . 
In the range 10-25 GHz the acceptance ratio increases 



to 80 %. This is well above the expected 1-a value of 
approximately 68 %. 



C. Analysis 

The results displayed in Fig. Ql|] are puzzling. The el- 
ements of S(f) are correlated over a frequency range T. 
They do not follow a Gaussian distribution. There are 
non-trivial higher-order correlations. On the other hand, 
there exists no discernible correlation among the rescaled 
Fourier coefficients Zk, and these are consistent with a 
bivariate Gaussian distribution for Sk/ y (xk) ■ How is it 
possible that a non-Gaussian distribution becomes Gaus- 
sian after Fourier transformation and rescaling? 

We first show that rescaling removes binary correla- 
tions. For simplicity, we do so for the case of a contin- 
uous frequency / ranging over the entire real axis. For 
clarity, we distinguish the ensemble average (indicated by 
an overbar) from the running average over the spectrum 
of a single realization (indicated by angular brackets). 

Without rescaling, the distribution of S{k) would ob- 
viously not be Gaussian. Moreover, rescaling does in- 
deed remove all correlations between pairs of S-matrix 
elements. To see this, we calculate the correlation func- 
tion of two Fourier-transformed S'-matrix elements, us- 
ing the translational invariance of the two-point correla- 
tion function [S fl (/i)5 fl *(/ 2 ) = S fl (/i + x)S ti *(f 2 + x) = 
.92 (/1 — f 2 ) for all real x]. That gives 



S(k 1 )S'(k 2 ) = 2n6(k 1 -k 2 ) sa(fci) 



(39) 



The Fourier transform 52 > of g 2 determines only the 
average value of |S(£;)| 2 ; pairs of Fourier-transformed S'- 
matrix elements with different arguments are uncorre- 
cted. We Fourier-transform S(k)/ \f]h(W) back to the 
frequency domain and find that the correlation function 
of a pair of Fourier-back-transforms is a delta function 
in frequency. Thus, the binary correlation has been re- 
moved by rescaling. Put differently, we may consider 
the quantities S{k)/ yjg 2 (k) as Fourier transforms of S'- 
matrix elements that are pairwise uncorrelated. 

Correlations of higher order (involving more than two 
S-matrix elements) imply correlations of higher order of 
the elements of S fl (fc) and of S(k)/ \fg2{k). Such corre- 
lations are not removed by rescaling. But they are made 
irrelevant by the way in which the distribution of the Zk 
is sampled. That is done by considering the index k as 
a label only. Any relation to the time scale originally 
inherent in the Fourier transformation is lost. We simply 
order the Zk by size, asking how many occur in each size 
interval. That yields the distribution in Fig. [19J The de- 
pendence on k is scrambled. It is not possible from that 
distribution to reconstruct correlations that may have 
existed among its elements. 

These arguments do not explain why the distribution 
of S fl (^)/ \J (xk) is Gaussian. (For that we must resort 
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to the law of large numbers) . But they show why corre- 
lations that are known to exist among the S R (f) do not 
prevent a Gau ssian to emerge for the distribution of the 
S & (k)/^k). 



D. Goodness-of-Fit Test 

To test the quality of the fit of the theoretical autocor- 
relation functions in Eqs. (fl~4| and Eq. (1241 to the data, 
we developed a goodness-of-fit test. The test applies to 
uncorrelated data with an exponential distribution. As 
shown in Sec. II V Cl that condition is met by the rescaled 
experimental Fourier coefficients Xk/(xk)- We recall that 
the mean values (xk) are determined by fitting a small 
number of parameters. This renders the decision non- 
trivial whether the fit is compatible with the data. As a 
measure for the goodness of the fit we used the expression 



A/-1 



fc=0 



In 



Xk 
(Xk) 



1 



(40) 



which quantifies the difference between the M data 
points Xk and the best- fit value (xk) for the theoretical 
expressions. The quantity / is non- negative and vanishes 
exactly if Xk — (xk) for all k. The expression Eq. (|4TJ)) 
is a generalization of the x 2 test used for Gaussian data, 
see Chaps. 14 and 16 of Ref. [62j| . If Xk/(xk) has an ex- 
ponential distribution then the distribution of / is given 
by 



1 f°° 
P(I) = — da, 

2 ^ J~oo 



,ia(I+M) 



T(l+ia) 
(1 + l a) 



(41) 



For our test we approximated P(I) by a chi-squared dis- 
tribution x ^ with M degrees of freedom 



(M/2) M / 2 
T{M/2)I 



exp 



MI 
~2~T 



(42) 



Here, / is the expectation value of / and is given by 



/ = M7. 



(43) 



where 7 = 0.577216 is Euler's constant. The agreement 
between P(I) and ^( M ) is better than 2 %. The same 
measure / was used for a goodness-of-fit test in Ref. [l2| ■ 
The test procedure is illustrated in Fig. [20l Upon the 
definition of a certain threshold K (in the Figure: K = 
0.9) on the cumulative distribution function a limit R 2 K is 
obtained (in the Figure: R 2 K = 1.064) which may not be 
exceeded by /// if the GOF test is to accept the model. 
The limit is defined by 



AIx {M) {I\I) = K . 



(44) 



The value of K quantifies the confidence into the test in 
the sense that 1 — K is the probability to make a wrong 
decision by rejecting a valid theory. We choose K = 0.9. 




FIG. 20: The probability distribution function (PDF, upper 
panel) and the cumulative distribution function (CDF, lower 
panel) of the distance measure I /I for M — 800. The solid 
lines in both panels correspond to the analytic model Eq. (|41|l , 
the results of a Monte Carlo simulation are shown in the upper 
(lower) panel as a histogram (circles). For a threshold of K = 
0.9 on the CDF a limit (horizontal dashed line) of B? K = 1.064 
is imposed on I /I (vertical dashed line) which may not be 
exceeded if the GOF test is to accept the model. 



In the fitting procedure the Fourier transforms of 
the theoretical autocorrelation functions of the GOE 
(Eq. (O), 01 the GUE (Eq. ([19])), and for the case of par- 
tial T-invariance violation (Eq. (I26|) h respectively, were 
fitted to the experimental Fourier coefficients. One ex- 
ample of such a fit is shown in Fig. 1211 In the upper four 
panels the results are compared to the data (dots) in the 
frequency domain. The lower panel shows the experimen- 
tal Fourier coefficients (dots) together with the best fits 
of the GOE (solid) and the partial TRIV (dash-dotted) 
result. We observe that the experimental Fourier coeffi- 
cients scatter widely around their average; the scatter is 
too large to directly arrive at a conclusive decision using 
the GOF test described above. Therefore data were taken 
from a total of six different realizations of the experi- 
ment. Moreover, the theoretical curves lie very close to 
each other. This corroborates the need of a goodness-of- 
fit test. This GOF test in conjunction with the enlarged 
data basis leads to the decision that within the confidence 
threshold i^T = 0.9 only the expression for partial TRIV 
describes the data. The GOE and GUE expressions are 
ruled out. 



V. SUMMARY 

We have measured reflection and transmission ampli- 
tudes in chaotic microwave billiards with two antennas. 
The measurements were performed in the regimes of iso- 
lated and of weakly overlapping resonances. Both a T- 
invariant system and a system with partially violated 
Tinvariance were investigated. The latter was realized 
by placing a magnetized ferrite within the microwave bil- 
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FIG. 21: Top: Comparison of autocorrelation functions in 
the frequency domain for a single realization in the range 
24-25 GHz, i.e. T/d « 1.14, and for B = 190 mT. The 
four panels display the results for Cn, C12, C21, C*22. The 
discrepancy between data (dots) and VWZ (solid) or complete 
TRIV (dashed) is more pronounced than for the model for 
partial TRIV (dash-dotted) with £ = 0.202. Bottom: C12 in 
the time domain (same key, but for clarity without the result 
for complete TRIV). 



Hard. The measurements yielded the moduli and phases 
of all four elements of the scattering matrix S(/) in a 
range of frequencies / limited by the requirement that 
only one vertical mode be excited in the billiard. The fre- 
quency range was divided into intervals of 1 GHz width. 
Within each interval, statistical measures for S-matrix 
fluctuations like the Fourier transform of the S'-matrix 
autocorrelation function, the distribution of S'-matrix el- 
ements, or the elastic enhancement factor were deter- 
mined from the data. 

We compared the results with theoretical expressions 
based on random-matrix theory. For T-invariant systems 
and for systems with full violation of Tinvariance these 
were given in Refs. [26|, [3(| [HJ while for systems with 
partial T-invariance violation they had to be calculated. 
This was done by extending the existing supersymmetry 
approach. The parameters of the theory are the trans- 
mission coefficients T\ and T2 for the two antennas, the 
parameter r a b s describing absorption in the billiard, and 
the parameter £ for the strength of T-invariance viola- 
tion. Starting values for these were partly obtained di- 
rectly from the data, but final values were always deter- 
mined from fits of the RMT expressions to some of the 
experimental measures. 



The large data sets taken made it possible to test 
the theoretical expressions with unprecedented accuracy. 
The outcome of these tests is recapitulated in Table U In 
particular, we used the following stringent tests, 
(i) Goodness-of-fit (GOF) test. For the T-invariant sys- 
tem, that test accepted the fit of the theoretical result for 
the Fourier-transformed autocorrelation function to the 
experimental data in all 1 GHz intervals. In the case of 
T-invariance violation, the GOF test was applied to data 
fits of RMT expressions for all three cases, i.e., the one 
for T-invariant systems, the one for systems with com- 
plete and the one for systems with partial T-invariance 
violation. The results of these fits are summarized in Ta- 
ble |U The fit is accepted in just 7 frequency intervals for 
the first case, it is rejected in all but two intervals for 
the second case, and it is accepted in all but one interval 
for the third case. In each 1 GHz window the fitted val- 
ues of £ agree with the values determined in Sec. IIIIBI 
We conclude that the GOF test is a powerful tool to un- 
cover the small effects of partial T-invariance violation 
on S-matrix fluctuation properties. 

TABLE I: Results of the GOF test for the billiard with ferrite 
at B = 190 mT. In each case the first row gives the lower 
boundary of the 1 GHz frequency interval used for the anal- 
ysis. The second (third) row indicates whether the GOF test 
for the autocorrelation function for Tinvariance (full violation 
of Tinvariance, respectively) was accepted. This is indicated 
by a bullet. If both expressions are accepted, no conclusions 
can be drawn and the column is marked by "o" signs. Re- 
sults rejected by the fit are indicated by "— ". The fourth row 
shows similarly acceptance or rejection of the expression for 
partial TRIV. 
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(ii) We inserted the fitted parameters into theoretical ex- 
pressions for the distribution of the diagonal S-matrix 
elements and, in the case of T violation, for the elastic 
enhancement factor. The results agreed well with the 
data. To the best of our knowledge, this is the first time 
that the elastic enhancement factor as function of the pa- 
rameter £ has been investigated in such detail over such 
a large frequency range. 

As an additional test, we extended our measurements 
beyond the frequency range where only one vertical elec- 
tric mode in the resonator is excited. When two such 
modes are excited, the resonator does not simulate a 
quantum billiard. The fit of the theoretical result based 
on random-matrix theory to the experimental data is re- 



17 



jected by the GOF test. That was shown in Ref. [12j 
and is not reproduced here. We have numerically sim- 
ulated the S-matrix correlations for that case under the 
assumption that the two modes of vertical excitation do 
not interact, using for the random Hamiltonian ensemble 
sets of real symmetric matrices consisting of two diago- 
nal blocks. Such matrices are commonly used to mimic 
spectral properties of "/"-invariant chaotic systems with 
some underlying symmetry. It was shown in Ref. [l2j 
that our simulations qualitatively reproduce the experi- 
mental autocorrelation function. The failure of the GOF 
test shows that our testing procedure is sensitive to the 
existence of such symmetries. 

We conclude that the theoretical expressions for the 
S-matrix correlation functions, for the distribution of S- 
matrix elements, and for the elastic enhancement factor 
based on random-matrix theory, are in excellent agree- 
ment with data measured on chaotic microwave billiards, 
both for the T-invariant case and for the case with partial 
7~-invariance violation. Our work constitutes the most 
stringent test of the statistical theory of quantum chaotic 
scattering yet done. The success in the case of partial T- 
invariance violation shows that the strength parameter 
£ can be determined reliably from scattering data. This 
is important in cases where that parameter cannot be 
reliably obtained theoretically from a dynamical calcula- 
tion like, for instance, the semiclassical approximation. 
The largest achieved values for the T-invariance viola- 
tion strength parameter £ equals 0.3. Numerical calcu- 
lations show that for this value the spectral fluctuations 
of the Hamiltonian H for the closed resonator defined in 
Eq. {TT} almost coincide with those of the GUE [Hj]. We 
also found that for £ = 0.4 the y do not differ significantly 
from those presented in Ref. [33], where the conclusion 
was drawn, that complete T breaking is achieved. How- 
ever, even for £ = 0.4 the value of C CIOSS (0) is still far 
from zero. This shows that C croS s(0) is a particularly 
suitable measure of the strength £ of T violation. 
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over the Grassmann variables which involves determining 
terms of highest (in our case eighth) order in the anti- 
commuting variables, is rather difficult. In Ref. [26| the 
integration was done with Efetov's original parametriza- 
tion for the i nteg ration measure and a method first devel- 
oped in Ref. (64|. The starting point was the generating 
functional 



Z(e) = Detg" 1 D + J(s 



exp ( — Trg In 



D + J(e) 



(45) 



Here, Detg and Trg denote the graded determinant and 
trace as defined in the supersymmetry formalism of 
Ref. [30| . The inverse propagator D and the matrix J 
are 47V x AN matrices. Introducing the matrix 



we have for the inverse propagator D 

1 



D = [f - H)I 4 + -eL + iQL 



(46) 



(47) 



Here ft = £ c ft c and L PP ' = (-l) p+1 8 pp >8 aa ' with 
p,p' = 1,2 and a, a' = 0, 1 is the diagonal supermatrix 
that distinguishes between the advanced (p = 1) and re- 
tarded (p — 2) parts of D. The index a — denotes the 
commuting, a — 1 the anticommuting components. The 
quantity e equals the difference of the arguments of the 
S'-matrix elements S ab (f - e/2) and S* b (f + e/2). The 
matrix J(e) is given as 

2 

■V({4&M&}) = ^EE^^abWV (48) 

a,b j = l 

where the matrix with entries 

I$(j) = (-l) 1+a S pp ,8 aa 'S„ (49) 

is the projector onto the p — j block. With Eq. © we 
obtain 
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Appendix 

For the derivation of the analytic expressions for the 
autocorrelation function and the cross correlation coef- 
ficient in Eqs. ([24]), ([26]) and ([27]), we used Efetov's su- 
persymmetry approach in the form of Ref. [26| . Due to 
the symmetry-breaking term in Eq. (|11[) . the integration 



S ab (f-e/2)S* ab (f + e/2) 

= 4Trft a Z3- 1 (/-e/2)ft fc [Z3- 1 (/ + e/2)] t , (50) 

where Tr denotes the trace over the index fi of the res- 
onator modes, and it can be checked that 



S ab (f - e/2)S* b (f + e/2) = 



U 2 



rlF 1 F)f 2 
Ut ab Ut ba 



Z(s) 



(51) 



e=0 



Averaging over the ensemble H becomes feasible when 
we write the generating functional as a Gaussian super- 
integral 



Z(e) = j P$ exp I l - ' 



^ pr ? 



(£> + J(e)) $ 



(52) 
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over an eight-dimensional supervector $. The matrices 
D and J have been extended to 8iV x 8N supermatrices, 



as 



and 



D 



3 = J{{e$ 1 }, 




(S)2, 



+^({4t )1 } ! Rr}) 



(A)2-, 



10 
l) 

' 1 
-1 



(53) 



(54) 



where e^f and e^"' are the parts of that are sym- 
metric and antisymmetric in the indices a and b, respec- 



tively. The indices r,r' = 1, 2 in Eq. (|52|) arise due to the 
doubling of dimension. The eight-component supervector 
<!> is given in terms of the four-component supervector <pp 

and its adjoint = <fls, where = (-1) (q+1)(1+p) <5^' 



<J> = (j)T s 



(55) 



The symmetrized form of the autocorrelation function, 
\ (S ab (f - e/2)S a * 6 (/ + 6/2) + S ba (f - e/2)S* ba (f + e/2)), 
is obtained by choosing the plus sign, that of the unnor- 
malized cross-correlation coefficient 9it({S a b(f)S ba (f))) 
by setting e = and choosing the minus sign in 



d 2 



4 de {S)1 de (S)2 

Ut ab Ut ba 



Z( £ ( s \0) 
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4 5e (A)1 5e (A)2 

Ut ab Ut ba 
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(56) 



e=0 



Ensemble averaging over the Hamiltonian H of the gen- 
erating functional Eq. (|52|) yields 



J 



<Z(e)> = / 2?$e l£(S) exp 




fI+-eL + iQL + 3(s) $ 



(57) 
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where 
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(58) 



with the supermatrix S^y r , = J^^^prM^p'r'it 1 ) an d 
T 3 r / = (— l) r+1 S rr >. The matrix L is eight-dimensional. 
It results from the doubling of dimension, L rr > = L5 rr ' ■ 
The quartic dependence of C(S) on $ is eliminated with 
help of a Hubbard-Stratonovich transformation. After 
expanding the resulting exponent in the large N limit 
in the small quantities e and —S- the remaining integral 
reads 



Here, trg denotes the supertrace over the indices p, a, r, 
and Trg includes the summation over the level index \x. 
The 8x8 supermatrix Q contains the commuting and the 
Grassmannian integration variables. For the integration 
Efetov's original parametrization of the Q matrix was 
used. 



(Z(e)} = / VQe 



JC eS (Q) e iC. IC (Q;J) 



(59) 



where £ e s(Q) = C tlcc (Q) + C ch (Q) and with X c = irNv 2 
(c.f. Eq. 0321) 
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i£ STC = --Trg In 




Q 



3(e) . (60) 



Before performing the integration (Z(e)) needs to 
be differentiated twice with respect to b and e ba at 
e[ b = and e{ a = for j = 1,2, c.f. Eq. (|55]l. 
Only terms of second order in e J ab and £ J feQ survive. 
Since J(e) is linear in e a b, these are obtained from the 
second order term in J of exp ^— -Trg In (l + M J 



which equals 



with M = Sr 
Eq. (HID) 



(-|Trg (m3\ ) + ±Trg (MJMj) 
^1 + . Then we finally obtain (see 
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4A 

(IT (-l) j ) 



(l + ^LQ 



h(l) trg 



All. 
A 



h(2) S ab (l±l) + a^b 



(61) 



For the integration over the Grassmann variables, we pro- 
ceeded as in appendix B of Ref. (2(|. The result is a 
threefold integral. For arbitrary values of the transmis- 
sion coefficients T c , the expression for the autocorrelation 
function for the case a ^ b is obtained from that given 
in Eq. (2) of Ref. [5l| by including in the integrand the 
factor exp (i^j (A1A2 — Ao)) arising for non-vanishing e 
from the last term in the first line of Eq. Here, 



Ao,Ai, A2 are the integration variables. It is straightfor- 
ward to compute the autocorrelation function for the case 
a = b by proceeding as in appendix B of [26j . The result 
for the cross-correlation coefficient is obtained by multi- 
plying the second and the third rectangular bracket in 
Eq. (2) of [5l| by (—1). All this yields for the autocor- 
relation function and the cross-correlation coefficient the 
expressions given in Eqs. (IM1) . (f2l)|) and (|2"7|) . 
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